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I.  A33TRACT 


‘neoretioal  efforts  at  understand ing  the  oand  structure 
of,  and  hence  electron  transport  in,  alloy  se.niconductor 
.naterials  are  descrioed.  In  particular,  calculations  of  oand 
structure  and  transient  and  equilibrated  velocity-field 
relationsnips  are  being  performed  for  Gj^n(T-^A3  and 
Qatftf-«7V(f-y>  naterials  which  nay  oe  suitaole  for  fast 
logic  or  nicrowave  devices.  Methods  for  describing 
pseudopotential  variations  and  chemical  disorder  in  alloys 
are  discussed . 


II.  INTRODUCriGM 

1 7 

Approxinate  calculations  of  velocity-field  curves  have 
indicated  tnat  the  alloy  se.nicond uc tor s Ga  ^In  -^As  and 
Ga  yyln  53AS  g,P  4.  which  can  be  epitaxially  natched  to  GaAs 
or  InP  substrates,  have  a higner  low-field  nooility  and  peak 
velocity  than  GaAs.  These  properties  nay  lead  to  i.nproved 
nicrowave  or  fast  logic  devices,  altnougn  experimental 
results  on  GalnAsP  microwave  FET's  nave  been 
disappointing  . 

Tne  work  covered  oy  this  grant  is  an  attempt  to  better 


' 1 
i 


describe  the  transport  properties  of  these  alloy  materials. 
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To  improve  tnis  description,  it  is  first  necessary  to 
examine  tneir  pand  structures  in  greater  detail  tnan  nas 
been  done  Pefore,  Consideraole  difficulty  arises  in  sucn  an 
invest! gation , However,  and  some  of  tne  difficulties  are 
descrioed  in  tnis  '•eport. 


111.  Band  Structures  for  IIT-V  Semiconductors 


A conputer  progran  was  written  and  successfully  applied 
to  tne  calculation  of  oand  structures  for  III-V  compounds 
using  nonlocal  pseudopotentials.  This  initial  pnase  of  the 
project  was  undertaken  as  a test  of  conputer  algoritnns, 
with  tne  intent  of  reproducing  tne  results  of  Cdelikowsxi 
and  Conen1.  Tne  oand  structure  calculation  we  performed 
generates  two  types  of  information: 

1.  Band  energies  as  functions  of  wave  vector.  From 
tnis  information,  one  can  extract  tne  following 
material  properties  needed  for  transport  calculations : 
A.  F undamental  gap , useful  for 

cnarac ter i zation  of  heteroj unctions  and 
rougn  estimates  of  ore akdcwn  field. 

3.  Band  shape  near  conduction  pand  ninima, 
needed  for  equations  of  motion  in  transport 
si nul ations . 

C.  Densi ty  of  s tates , needed  to  calculate 
scattering  rates. 


0.  Sub-band  gaps,  needed  for  simulating 
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intervailey  transfer. 


2.  Pseudo-wave  functions.  These  functions  are  used  in 
subsequent  calculations  to  yield  nany  necessary 
quantities,  including: 

A.  Intervalley  scattering  rates  for 
transport  simulations. 

B.  Coupling  strengths,  which  could  be 
useful  in  tne  developnent  of  a conplete 
theory  of  alloy  scattering. 

C.  Matrix  elements  for  a wide  variety  of 
other  transitions,  such  as  pnotoaosorpcion 
and  impact  ionization. 

^e  verified  our  computer  progran  by  calculating  the  band 
structure  of  GaP  at  symnetry  points  r,  L and  , and 


comparing  with  published  results: 
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17.  Band  Structures  for  Alloy  Semiconductors 


alculations  of  alloy  oand  structures  by  otne 


nave 


involved  several  approximations 


For  exanple,  all  have  used 


local  pseudopotentials,  wnicn  results  in  equal  treatment  of 


the  wave  function  components  of  various  angular  momenta 


energy  dependence  of  the  pseudopotential  operator 


Also  tn 


These  calculations  appear  to  be 


nas  been  ignored 


inconsistent  with  experimental  data 


indicating  that  a 


nonlocal  operator  may  be  necessary  for  alloys 


Tnese 


inconsistencies  will  be  discussed  later.  Reinforcement  fo 


this  statement  is  furnished  by  the  report 


tnat  nonlocal 


contributions  are  important  in  the  calculation  of 


intervalley  phonon  coupling  constants 


even 


compound  semiconductors 


In  an  "empirical"  band  structure  calculation,  i.e 


one 


in  which  tne  pseudopotential  parameters  are  adjusted  to  give 


reasonable  agreement  between  calculated  and  measured 


spectral  properties 


the  nonlocal  netnod  requires  ncre 


local  nethod 


fortunately 


parameters  than  tn 


tnese 
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nonlocal  parameters  nave  already  oeen  determined1  for  a 
large  nunoer  of  III-V  compounds  in  a nanner  tnat  is 
apparently  sufficiently  accurate  for  tne  deter nination  of 
input  paraneters  for  transport  simulations . 

As  nave  developed  techniques  for  interpolation  of  the 
local  par t of  the  pseudopotential  across  alloy  composition 
ranges.  Tne  empirically  determined  form  factors  V(k)  for 
eacn  compound  at  its  equilibrium  lattice  constant  are 

p 

subjected  to  a nonlinear  tr ansfor nation  to  remove  tne 
screening,  then  fitted  to  the  tnree-parameter  model  function 

V(k)=AU-3)sxp(-C(k-t3) ) . 

This  model  function  is  then  evaluated  at  wave  numaers 
appropriate  to  the  lattice  constant  of  tne  alloy,  and 
renormalized  with  respect  to  tne  atomic  volume  of  tne  alloy. 
The  forn  factors  for  the  two  or  four  compounds  (for  a 
ternary  or  quaternary  alloy,  r espectively)  are  linearly 
interpolated  according  to  tne  alloy  composition  ratio(s), 
and  screening  is  reapplied. 

Two  techniques  were  developed  for  interpolating 
the  nonlocal  part  of  the  pseudopotential  in  alloy  systems. 
Unfortunately , neither  nas  yet  given  results  wnicn  can  be 
accepted  with  confidence. 

In  the  first  tecnnique,  tne  cuoes  of  the  nonlocal  radii 
(3^)  were  linearly  interpolated;  tne  st-engtn  parameters  a 
and  8 were  then  renormalized  wi  tn  to  H^,  and 
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subsequently  interpolated  . For  a ternary  alloy, 

(8>alloysx(S)tia/t<alloy  * < 1 “x) (8)ftp/rtalioy * 
wnere  x is  the  composition  ratio. 

Tin i 3 netnod  is  inspired  physically  oy  one  fact  tnac  in  most 

of  tne  natrix  elener.ts  of  tne  nonlocal  operator,  tna 

nonlocal  parameters  enter  as  (a+*g)R^,  so  the  net  nonlocal 

strength  should  be  approximately  preserved  if  a and  6 are 

scaled  along  with 

Calculations  on  Ga  ^ln  -^As  using  tnis  scheme  resulted 

spin  orbit 


in 


bandgap  Ejs0.97  eV  (corrected  for 


splitting),  higher  than  the  measured  value  of  J.dO  eV  . 
Part  of  this  difference  nay  oe  attributable  to  neglect  of 
alloy  conposi tional  disorder  ir.  the  calculation,  which  uses 
the  "virtual  crystal"  approximation  (JC A).  A disorder 
correction  tern  can  oe  added  to  the  VC A energies,  which 

p 

results  in  a decrease  of  tne  calculated  gap  to  aocuc 
0.95eV.  As  conclude  that  either  our  VGA  calculations  are  in 
error,  or  the  published  results  for  the  disorder  effects  are 
inaccurate,  possibly  due  to  tne  fact  that  they  were  found 
from  a local  pseudopotential . 

In  order  to  verify  the  correctness  of  our  VC  A 
calcul ations , a second  nonlocal  pseudopotenti al 
interpolation  scheme,  which  involves  calculating  each  natrix 
elenent  separately  for  eacn  conpound  then  perforning  tne 
linear  interpolation  or.  the  nat-ix  ele.ier.ts,  was  tried: 

^’■V^lloy  = 
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x<k-3f  !H"'1*Sic-34>  ♦ (1-x)<k-Gj  !i<-G  .> 

wnere  !k-G>  is  a plane  wave  basis  state. 

This  fornula  correctly  scales  the  nonlocal  terns  only  for 
tne  alloy  ato.nic  vclune,  not  for  the  change  in  screening. 
However,  screening  may  not  be  inportant  for  the  nonlocal 
part  of  tne  potential1. 

Renor nal ization  of  R to  account  for  the  lattice  constant 

change  has  been  considered,  but  was  not  done  oecause  tne 

available  netnods  (sucn  as  a quadratic  variation3  or  tne 

linear  "touching  spheres"  criterion1)  would  require  previous 

Knowledge  of  the  band  structure  for  each  alloy  composition. 

1 A 

Since  tne  results  are  relatively  insensitive  to  the  value 
of  R as  long  as  a and  8 are  chosen  appropriately,  R was  left 
unchanged  fron  its  value  in  the  binary  compound.  Tne 
results  for  GalnAs  using  this  second  interpolation  paradigm 
were  nearly  identical  to  the  previous  results. 

On  the  oasis  of  the  above  results,  it  was  decided  that 
new  calculations  of  the  disorder  effect  should  be  performed, 
using  the  nonlocal  pseudopotential,  in  order  to 

1.  satisfy  the  goal  of  calculating  accurate  suo- 
band  gaps  for  alloys,  and 

2.  deternine  whether  the  alloy  pseudopotential 

is  oeing  properly  formulated,  thereoy 

establishing  a confidence  level  for  tne  matrix 
elenents  that  will  be  derived  from  tne  wave 
functions,  even  though  there  are  presently  no 
plans  to  correct  the  wave  functions  for  disorder 


pa  at  i 

effects.  Until  tne  energy  band  results 


(VCA+disorder) 


are  nade  to  agree  witn 


neasur ernents , there  is  no  guarantee  that  tne  VCA 
pseudcwave  functions  are  accurate. 


Tnere  is  considerable  confusion  in  the  literature 

•5  10 

regarding  the  disorder  effect.  Some  authors^’  appear  to 
have  obtained  good  agreenent  with  experiment  Oy  ignoring 
disorder,  wnile  otners11'1^  nave  stressed  the  importance  of 
disorder,  particularly  for  anionic  nixtures  sucn  as  GaAsP. 


A local  pseudopotential  calculation 


of  the  disorder- 


induced  sowing  of  the  direct  gap  in  GaAsP  shows  it  to  oe 

nearly  equal  to  the  observed  bowing,  inplying  that  tne  VCA 

ocwing  should  be  nearly  zero.  However,  our  local 

pseudopotenti al  calculation  of  VCA  oowing  gives  -J.46  eV. 

It  is  commonly  agreed  that  disorder  is  relatively 

unimportant  in  AlGaAs  and  AlGaSd  J,  out  there  is  seme 

debate  regarding  GalnAsP.  Some  treatments1*^  assune  that 

disorder  affects  tne  entire  Srillcuir.  zone  with  tne  sane 

strength  while  others11  claim  that  some  points  are  ncre 

1 4 

strongly  affected  than  others.  Altarelli  states  that 
disorder  is  mere  important  for  the  conduction  band  than  for 
the  valence  band,  while  Saldereschi  and  Mascnke13  find  that 
the  reverse  is  true.  Tne  latter  wo r<  also  conflicts  with 
Stroud's  clain11  that  disorder  always  reduces  the  gap 
relative  to  tne  VCA  value.  Tne  evaluation  of  oand 
structure  calculation  techniques  for  alleys  is  hindered  oy 
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tne  prolif jration  of  incorrect  pseudopotential  parameters  1 ^ 
and  oy  tne  spread  in  experimental  data. 


-•  Jertsi cy  c f States  Calculations 


To  reduce  conputation  tine  during  transport  simulations, 
a two-band  secular  equation  nodel  was  applied  to  tne 
"camel's  oac*"  conduction  oand  structure  at  tne  X symmetry 
point,  and  expressions  ware  derived  for  energy  dispersion 
relations  and  density  of  states.  Tne  nodel  is  used  ay 
fitting  to  several  points  frcn  the  full  pseudcpotential 
calculation. 

For  GalnAsP  composi tions  rich  in  As  and  In,  tne  A 
•nininun  is  expected  to  be  at  X and  the  dispersion  relation 
is  nearly  paraoolic.  For  conpositions  rich  in  Ga  and  P,  tne 
mininun  shifts  toward  f 3nd  tne  oar.d  shape  is  approximately 


EU)sA.<2+Bi<2-[  ( A/2)2«-AQA.<2]  1 /2 
where  and  «t  3re  the  components  of  a (measured  with 
respect  to  X)  parallel  and  perpendicular , r espectively , 
to  the  line.  Aaft^/Sm^,  3=l72/<?:nt. 

The  density  of  states  is  given  by 


ra. 


"2  *2/?a" 


( v'x«-y  - a.,(x-y)/^y) 


where  xs  Ao-2(£  -E), 

y/rf5S2*4i  u-t;) 


u_j  = unit  step  function. 

Tne  paraneters  A,  AQ  , £c , A and  B are  cotained  oy  fitting 
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to  calculated  band  snapes. 

V.  CJMCLJSUdS 

A ccnputer  progran  was  developed  for  the  purpose  of 
calculating  band  structures  for  III-/  compounds  using 
nonlocal  pseudopo tenti al s . Upon  verification,  tnese 

tecnniques  were  extended  to  band  structure  calculations  of 
III-'/  alloys,  but  the  results  cannot  yet  oe  reconciled  witn 
experinental  data.  Additional  calculations  are  now  being 
perforned  for  tnat  purpose. 

A large  body  of  expari nent3l  measur enents  is  availaole 

for  ccnparison  with  calculated  results,  including  energy 
21  Id 

gaps  , effective  nasses  , and  high-field  electron 
transport  data  . 

A possible  cause  of  disagreement  between  tn.ecry  and 

experinent  is  the  preferential  clustering  (also  called 

segregation  O’*  snort-range  order)  which  nas  oeen  ooserved  in 

alloys.  Tnere  have  been  theoretical  attenpts  aimed  ac 

o p 

predicting  the  extent  of  this  pnencnenon  , and  segregation 

has  Oeen  ooserved,  for  example,  in  Hr  id g man-grown  ouIk 

GalnP^.  The  influence  of  short  range  order  on  oand 

2 4 

structures  has  oeen  studied  for  binary  systems,  wnere  tne 
effect  was  described  as  "striking,"  and  tne  extension  to 
pseudo-oinary  systems  like  GalnAsP  should  be  feasible.  The 


effect  of  segregation  as  an  additional  scattering  mechanism 


(pernaps  to  be  used  in  place  of  tne  usual  "alley  scattering" 
tern  which  was  derived  cn  tne  assunption  of  no  segregation) 
nas  oeen  introduced  into  electron  transport  simulations^. 

*/e  consider  tne  following  topics  to  oe  worthy  of  further 
st Jdy : 

1.  Disorder  contr ioutions  to  the  band  structures 
of  alloys.  Tne  validity  lor  invalidity)  of  the 
virtual  crystal  app-ox imat ion  snould  be 
definitively  estaolished. 

2.  Alley  scattering.  Accurate  methods  snould  oe 
estaolisned  to  replace  the  educated  guesses^ 
which  nave  oeen  used  to  deternine  the  strength  of 
this  process. 

3.  Interpolation  methods  for  pseudepo tential 
parameters  snould  be  studied  more  thoroughly. 

Tne  various  ncdel  potentials,  such  as  square  well, 
Gaussian  well,  and  1/r,  should  oe  ccnpared 
systematically. 

4.  Short  range  order  in  alloys  should  oe 
investigated  experimentally  for  various  crystal 
growth  conditions,  and  the  results  snould  oe 
incorporated  into  Oand  structure  and  transport 


calculations . 
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SAMPLE  ALLOY  SEMICONDUCTOR  BANDSTRUCTURE  CALCULATION 
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Transient  velocity  characteristics  of  electrons  in  GaAs  with 
r-L-X  conduction  band  ordering*’ 

S.  Kratzar  and  Jaffray  Fray 

School  of  Electrical  Engmeenog.  Cornell  Uniremty.  Ithaca.  Sew  York  1 4853 
(Rse*ved  to  February  1971;  accepted  for  pubtocaiioo  I March  1976) 

The  transient  velocity  characteristics  of  electron*  in  OaAs  have  been  calculated  for  uniform  electric  field* 
el  room  temperature.  Three-level  bend-etructurc  data  have  been  uted  with  a Monte  Carlo  umulauon 
program  to  predict  electroo  drift  velocity,  temperature,  and  valley  population  fractions  as  functions  of 
position  and  electric  field  Results  lor  two  three-level  hand-structure  models  are  compared  with  those  lor 
the  previously  accepted  two-level  model.  Tranaient  and  steady-state  properties  are  seen  to  be  sensitive  to 
some  material  parameters  whose  values  arc  soil  in  dispute. 


FACS  cumbers:  72.S0.Ey.  72.20.Ht 

The  performance  of  very  small  electron  devices,  such 
as  microwave-field-effect  transistors,  may  well  be  af- 
fected by  the  transient  electron  transport  properties  of 
the  semiconductor  materials  of  which  they  are  made. 
This  paper  describes  an  investigation  of  some  of  the 
transient  properties  of  electrons  in  GaAs,  using  recently 
published  data1  concerning  the  three- level  band  struc- 
ture of  this  material  and  the  processes  that  scatter 
electrons  among  the  levels. 

I.  INTRODUCTION 

The  necessity  for  this  work  arises  from  new  informa- 
tion on  the  band  structure  of  GaAs.  Previous  calcula- 
tions of  the  transient  response  of  electrons  in  this  mate- 
rial1"* used  a two-level  band-structure  model  which,  in 
the  light  of  new  experimental  results,  is  now  taken  to  be 
incorrect.  This  two- level  model  (hereafter  called  model 
T)  placed  the  X minima  of  the  conduction  band  at  a lower 
energy  than  the  L minima,  and  the  L-X  separation  was 
assumed  to  be  large  enough  so  that  scattering  to  the  L 
valleys  did  not  occur.  In  contrast,  the  new  data,  ob- 
tained by  electroreflectance5  and  other  techniques/  have 
been  interpreted  recently  to  show  that  the  correct  se- 
quence, in  order  of  increasing  energies,  is  T-L-X; 
further,  the  intervalley  separations  are  small,  so  that 
scattering  among  all  valleys  must  be  considered. 

Several  authors1’*  have  reported  theoretical  calculations 
which  imply  that  the  new  band- structure  model  can  be 
consistent  with  the  observed  steady-state  transport  pro- 
perties of  GaAs.  We  recall  that,  as  early  as  1970,  Faw- 
cett, Boardman,  and  Swain'  expressed  consternation 
over  the  fact  that  some  experimental  results  were  in- 
consistent with  model  T. 

It  is  possible  to  study  the  transient  behavior  of  elec- 
trons in  small  devices  by  solving  the  differential  equa- 
tions describing  conservation  of  energy  and  momen- 
tum. However  solution  of  these  equations  becomes 
simple  only  when  some  extreme  simplifications,  9uch 
as  the  use  of  relaxation  times  derived  trom  steady- 
state  distribution  functions  are  made.  The  accuracy  of 
the  results  using  this  technique,  therefore,  is  question- 
able. The  preferred  method  is  the  more  accurate  but 
more  costly  one  of  Monte  Carlo  simulation,  described 
elsewhere. 7 The  transient  properties  are  calculated 


‘’Work  supported  by  the  U.S.  Air  Force  Office  of  Scientific 
Research. 


here,  using  this  method,  (or  electrons  in  GaAs  under 
the  Influence  of  a uniform  electric  field.  The  scattering 
processes  Included  in  the  transport  model  are  acoustic 
phonon  through  deformation  potential,  polar  optic  phonon, 
intervalley,  ionized  Impurity,  and  (for  the  L valleys 
only)  nonpolar  optic  phonon. 

II.  MATERIAL  PARAMETERS 

Values  of  the  electron-lattice  interaction  parameters 
of  GaAs  for  this  calculation  were  obtained  by  critically 
examining  all  available  data,  which  for  some  para- 
meters cover  a wide  range.  In  Table  1,  column  T re- 
presents the  two- valley  model  used  for  previous  trans- 
ient calculations.  Column  B shows  values  obtained 


TABLE  1.  Parameters  (or  two-level  IT'  and  three-level 
<B,C)  models  for  the  band  structure  of  GaAs. 


Param- 

eter 

Units 

T* 

C 

Alternatives 

ml,c 

”»o 

0.  322 

0.128 

0.35,  c 0.  IT1 

m*.c 

”»o 

0. 35 

0.58 

0.58 

0.  35, c 0.  33* 

mt.D 

w„ 

0.  222 

0.244 

0.  35,  c 0.  26d 

’”1.0 

m0 

0.35 

0.  58 

0.58 

0.35,  * 0.43d 

En 

eV 

0.  33 

0.30 

0.  38,  c 0.  285  • 0.  29 1 

Erx 

eV 

0.36 

0.  522 

0.48 

0.40,  0 0.33,* 

0.485,'  0.36" 

£,(D 

eV 

7.0 

7.0 

7.0 

Et(L) 

eV 

9.2 

2. 24 

7.  O' 

E,  CY) 

eV 

7.0 

9.27 

3.0 

7.0' 

a(D 

eV( 

0.576  0.610 

0.610 

all) 

eV*1 

0.461 

0.461 

a IX) 

eV*1 

0 

0. 204 

0.204 

&TI 

10*  eV/cm 

10 

10 

2. 8C 

Drx 

10*  eV/cm 

10 

10 

10 

a,'  5* 

Du 

10*  eV/cm 

10 

10 

1.8' 

DXx 

10*  eV/cm 

10 

7 

7 

10,  ' 1.  3* 

Due 

10*  eV/cm 

B^SS 

5 

5 

5. 5' 

‘Reference  7. 

‘Reference  i. 

•Reference  6. 

*1.  Topol,  H.  Neuman,  and  E.  Hess,  Phys.  Status  Solldl  B 
58,  K47  (1973). 

‘Reference  5. 

* T.  Sugeta,  A.  Mnjerfeld,  A.  K.  Saxena,  and  P.  N.  Robson 
(unpublished). 

*L.W.  James  and  J.L.  Moll,  Phys.  Rev.  183,  741  <19691. 
‘A.R.  Hutson,  A.  Jayaraman.  and  A.S.  Corlell,  Phvs.  Rev. 
155,  786  (1969). 


4064 


J.  Appi.  Phys.  46(7).  July  1678 


0021  8979/78/4907  4064801  10 


t£  1978  American  institute  of  Phytic* 


4064 


for  a three- level  model  by  Littlejohn  et  al. 1 in  a pro- 
cedure aimed  at  achieving  a “best  fit”  between  Monte 
Carlo  simulations  and  steady- state  drift  velocity  mea- 
surements. This  model  (B)  was  used  as  a starting  point 
in  our  investigation  of  transport  properties.  We  de- 
scribe below  the  reasons  for  considering  changes  to 
some  of  the  parameter  values  used  in  model  B. 


The  nonparabolicity  values  for  the  L and  X valleys  of 
model  B were  taken  to  be  isotropic  and  were  calculated 
from  pseudopotential  band- structure  predictions. 10 
Jacoboni  el  al.,*  in  their  Monte  Carlo  simulation  of 
high- field  transport  in  the  A valleys  of  Si,  also  used  an 
isotropic  nonparabolicity  constant,  but  they  obtained  the 
value  empirically  rather  than  from  a band- structure 


The  effective  masses  for  the  satellite  valleys  in  model  calculation.  However,  even  if  the  effective  masses  for 
B are  isotropic;  that  is,  the  conductivity  mass  m?  and  these  valIeys  were  isotroPic  at  the  symmetry  points, 

density-of-states  mass  mj  are  taken  as  equal.  However,  the  ^Parabolicity  constants  would  not  be  so.  The  band- 

electroreflectance  measurements5  and  pseudopotential  structure  curves,  by  convention,  cover  only  a small 
calculations10  indicate  that  the  L valleys  are  highly  an-  of  the  BriUouln  zone‘  For  the  L valleys,  only  long- 

isotropic.  Therefore,  in  model  C we  have  obtained  itudinal  variations  are  shown  (T  to  L,  but  not  L to  K or 

values  tor  ml  by  correcting11  the  data  of  Ref.  5 for  in-  L t0  ' The  nonparabolicity  value  obtained  from  these 

version  asymmetry  and  coupling  to  higher  conduction  data’  therefore>  is  appropriate  only  for  the  longitudinal 

bands.  component  of  k.  Paige,  “ in  his  analysis  of  the  L valleys 

of  Ge,  assumed  that  nonparabolicity  is  significant  only 

The  subband  naps  used  in  model  B were  taken  from  the  for  the  transverse  component  of  k.  James16  suggests 
most  recent  measurements  (£ri  by  high-pressure  tran-  using  a high-order  polynomial  in  kx,  kr,  and  kx  (rather 

sport,  El x by  electroreflectances5).  An  additional  than  Iftl)  to  account  for  the  anisotropy  of  the  nonpara* 

measurement,  of  Erx  by  optical  absorption, 13  was  not  bolicity.  This  topic  needs  further  investigation,  but  in 

used  in  the  formulation  of  model  B.  However,  this  ne-  the  absence  of  complete  calculations  we  have  retained 

glected  last  measurement  yields  data  subject  to  a higher  the  values  of  model  B. 
degree  of  experimental  control  than  that  obtained  by 

high-pressure  transport.  Further,  the  analysis  of  the  Finally,  we  note  that  Adams  el  al."  presented  a com- 

electroreflectance  structure  leads  to  an  uncertainty  of  plete  model  for  GaAs  material  parameters.  We  did  not 
30  meV  for  Erx. 11  Therefore,  we  have  used  what  we  use  this  model  because  it  includes  a value  for  w£ 

believe  are  correctly  weighted  averages  of  all  avail-  35'”o),  which  apparently  cannot  be  justified  by  any 

able  data  for  the  subband  gaps  in  model  C.  published  experiment  or  calculation.  In  the  trial-and- 

error  fitting  process  used  by  these  authors,  an  error 
The  acoustic  deformation  potentials  of  model  B for  the  in  one  parameter  can  lead  to  compensating  errors  in 
satellite  valleys  were  obtained  from  measurements  on  several  others,  so  that  nonphysical  parameter  values 
indirect-gap  materials.  Such  measurements  can  yield  may  result  even  while  the  end  results  of  the  calcula- 


a wide  range  of  values;  for  example,  those  for  the  X 
valleys  of  GaP  range14-16  from  3.7  to  55  eV.  Some 
authors1, • 18  have  attempted  to  calculate  the  acoustic  de- 
formation potential  from  derivatives  with  respect  to 
pressure  of  direct  energy  gaps,  taking  the  deformation 
potential  for  the  valence  band  as  zero.  However,  Rode19 
has  noted  that  this  procedure  may  be  very  inaccurate 
and  can  be  justified  only  because  this  scattering  mecha- 
nism is  relatively  unimportant  in  GaAs  at  room  tem- 
perature. In  the  interest  of  rigor,  therefore,  the  values 
of  acoustic  deformation  potential  for  the  L and  .Y  valleys 
in  column  C were  obtained  from  the  dielectric  bond 
calculations  of  Camphausen  et  al.,10  which  include  the 
strain  dependence  of  the  valence  band  edge. 

It  should  be  noted  that  the  strained- crystal  pseudo- 
potential method  of  Tsay  and  Bendow21  could  be  useful 
for  calculating  acoustic  deformation  potentials  of  GaAs. 
However,  these  authors  have  not  repcrted  some  of  the 
data  necessary  for  the  calculation,  viz,  deformation 
potentials  for  the  Indirect  gaps  and  the  ionization  poten- 
tial. Yee  and  Myers”  also  used  a nonperturbative 
pseudopotential  method  to  calculate  the  acoustic  de- 
formation potential  for  the  T valley.  The  computations 
are  relatively  simple  and  could  easily  be  extended  to 
the  satellite  valleys.  However,  the  additive  arbitrary 
constant  which  arises  in  their  band- structure  calcula- 
tion has  not  been  shown  to  be  volume  independent. 
Therefore,  as  Harrison”  points  out,  such  calculations 
could  be  used  to  find  deformation  potentials  of  energy 
gaps,  but  not  the  deformation  potential  of  the  conduction 
band. 


tions  may  agree  with  some  experimental  data. 

In  the  results  presented  below,  models  T and  B are 
compared  in  order  to  determine  whether  previous  GaAs 
device  analyses8,37  would  be  significantly  altered  by  the 
use  of  a three-level  material  model.  Models  B and  C 
are  compared  in  order  to  demonstrate  that  the  individ- 
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FIG.  1.  Fraction  of  total  electron  density  In  each  set  of 
satellite  valleys,  as  functions  of  electric  field,  for  models  T 
and  B.  .V»  10l£,  r»300  K.  For  model  C,  X valley  fraction  Is 
same  as  for  model  B;  L valley  fraction  is  nearly  equal  to  X 
valley  fraction  of  model  T. 
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FIG.  2.  Equilibrium  drift  velocity  vereua  electric  field  for 
models  T,  B,  and  C;  <V«  1011. 


ual  parameter  changes  described  above  have  a sub- 
stantial effect  on  the  predicted  transport  properties. 
However,  model  C shares  at  least  some  of  the  failings 
of  model  B since  the  interv  alley  coupling  constants  and 
mg  values  of  model  B,  which  were  chosen  to  give  rea- 
sonable steady- state  properties  in  conjunction  with  the 
other  parameters  of  this  model,  have  been  transferred 
to  model  C without  alteration. 

III.  STEADY-STATE  RESULTS 

To  obtain  some  information  on  the  differences  among 
models  T,  B,  and  C,  we  have  calculated  the  steady- 
state  transport  properties  using  each  model.  Figure  1 
shows  that  for  model  B the  X valley  population  is  quite 
small  for  fields  up  to  20  kV/cm.  However,  this  does 
not  mean  that  a two- valley  (r -L)  model  would  be  ade- 
quate for  low-field  Monte  Carlo  simulations,  because 
T—X—L  transfers  appear  to  be  important  in  models  B 
and  C.  Figure  2 shows  that  model  C has  a tower  peak 
velocity  and  higher  saturation  velocity  than  model  B, 
but  both  models  yield  results  that  fall  within  the  range  of 
experimental  data,  indicated  by  the  bars  in  Fig.  2,  up 
to  about  10  kV/cm.  At  higher  fields,  model  C gives  a 
velocity  that  diverges  from  the  results  of  model  B and 


Model  B 


FIG.  3.  Equilibrium  velocity  versus  electric  field  for  model 
B;  T-300  K,  .V-0,  10“,  and  10,T. 
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FIG.  4.  Electron  temperature  versus  electric  field  for  model 
B;  N » 101T.  Dashed  line  represents  an  empirical  formula  to 
be  used  for  FET  noise  analysis. 


from  experiment,  primarily  because  of  the  differences 
in  ml  between  the  two  models.  The  effect  of  doping  on 
the  steady- state  properties  of  model  B is  shown  in  Fig. 
3. 

Electron  temperature  as  a function  of  electric  field  is 
shown  in  Fig.  4,  for  model  B.  The  three-level  model 
(B)  yields  temperatures  somewhat  lower  than  does  the 
two-level  model  (T),  following  the  law  T,  = 300(l  + 1.8 
(£/5)2,4]  for  £<4.5  kV/cm.  Therefore,  GaAs  transis- 
tors should  have  lower  intrinsic  noise  figures  than  were 
previously  predicted  on  the  basis  of  the  two-level 
model. 22  In  fact,  the  latest  GaAs  FET  noise  results28  do 
show  this  property.  In  addition,  recent  experimental  re- 
sults for  GaAs  layers29  have  shown  a Tt — vs — electric- 
field  relationship  in  agreement  with  that  calculated 
here. 

IV.  TRANSIENT  TRANSPORT  BEHAVIOR 

Figure  5 illustrates  the  velocity  overshoot  phenome- 
non for  models  T and  B in  terms  of  velocity  versus 
distance  for  electrons  which  are  initially  in  thermal 
equilibrium  at  300  K.  The  field  strengths  of  10,  25,  and 
75  kV/cm  have  been  taken  as  covering  a range  of  values 
which  may  occur  in  field-effect  transistors  with  gate 
lengths  of  1 jam  or  less. 

Figure  6 shows  the  populations  of  the  satellite  valleys 
as  functions  of  distance,  for  models  T and  B.  At  the 


FIG.  5.  Instantaneous  velocity  versus  distance  for  models  T 
and  B;  101T,  £«  10,  and  25,  and  75  kV/cm. 
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FIG.  6.  Satellite  valley  fractions  versus  distance,  for  models 
T and  B;  AT-  iOn,  £-10,  25.  and  75  kV/cm. 


higher  field  strengths,  model  B Is  seen  to  exhibit  rapid 
transfer  to  the  L valleys,  via  V~L  and  r—.Y —I  transi- 
tions. After  the  L valley  population  has  become  signifi- 
cant, L—X  transfer  causes  the  X valleys  to  become 
populated. 

The  effects  of  doping  on  the  transient  properties  of 
model  B are  shown  in  Fig.  7.  As  previous  authors  have 
found  for  model  T,  impurity  scattering  has  very  little 


Model  B 


FIG.  7.  Instantaneous  velocity  versus  distance  tor  model  B; 
£-10  kV/om,  JV-0  and  10”. 
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FIG.  8.  Instantaneous  velocity  versus  distance  for  models  B 
andCj.V-lO1’,  £-10,  25.  and  75  kV/cm.  At  23  kV/cm,  the 
two  curves  are  Identical  to  within  computational  accuracy. 


effect  at  high  field  strengths  and  causes  a reduction  in 
peak  velocity  at  lower  field  strengths. 

Figure  8 shows  a comparison  of  transient  responses 
for  models  B and  C.  At  10  kV/cm,  the  peak  velocity  of 
model  B is  JO'I  higher  than  that  of  model  C,  but  at 
higher  fields,  the  peak  velocities  are  nearly  equal. 

V.  CONCLUSIONS 

The  transient  response  of  electrons  in  GaAs  has  been 
calculated  using  a now-discredited  two-level  band- 
structure  model  and  two  three-level  models.  A review 
of  available  data  on  material  parameters  has  revealed 
some  apparent  inadequacies  in  individual  parameters  of 
one  of  the  three-level  models,  which  had  been  optim- 
ized as  a set  of  parameters  to  yield  steady- state  pro- 
perties which  agree  with  experiment.  Steady- state  tran- 
sport simulations  suggest  that  the  other  three-level 
model,  which  was  arrived  at  by  consideration  of  the 
model  parameters  individually  rather  than  as  a com- 
plete set,  is  also  not  fully  satisfactory.  Nevertheless, 
the  transient  results  presented  herein  should  be  useful 
in  the  further  analysis  of  such  GaAs  devices  as  micro- 
wave-field-effect transistors. 

The  predicted  transient  properties  are  found  to  be 
quite  sensitive  to  changes  in  the  material  parameters, 
even  if  we  restrict  our  attention  to  those  models  that 
give  good  agreement  with  measured  steady- state  drift 
velocities.  While  this  fact  complicates  the  analysis  of 
tiny  electron  devices  when  several  material  models  are 
being  considered,  it  also  provides  an  opportunity  to 
select  from  among  these  models  when  experimental 
transient  results,  as  can  be  obtained  by  a microwave 
time-of-fllght  experiment, 50  become  available. 
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